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Screening of an external random potential by a two- 
dimensional (2D) Fermi liquid may be strongly dependent on 
the degree of its polarization. This is proposed as a possi- 
ble mechanism for the observed strong magnetoresistance of 
the 2D electron liquid in the magnetic field parallel to the 
electron plane. A Hartree-Fock calculation for the weakly 
disordered Hubbard model which qualitatively describes the 
experiments on the diagonal and the Hall resistivity in the 
finite-temperature metallic state in the high-mobility Si in- 
version layers is presented. 



The problem of a possible metal-insulator transition in 
two-dimensional (2D) electronic systems has attracted 
a considerable amount of attention recently M. The 
main reason for the excitement is that the existence of 
a metallic state at zero temperature (T = 0) would run 
against the intuition built on the weak-localization the- 
ory for non- interacting particles in 2D 0], and indicate 
that some new physics has resulted from strong electron 
correlations. The experiments are, however, unavoidably 
performed at finite temperature, and the issue whether 
the observed finite temperature metal would survive at 
arbitrarily low temperatures, or weak- localization would 
eventually take over, is presently hotly debated ||§1, [Q, 
Pi, 0. In either case, the finite temperature metallic 
phase shows some unusual transport properties and is 
quite interesting in itself. One such feature, which has 
often been invoked as a demarcation criterion for the 
metallic phase, is its strongly temperature dependent re- 
sistivity, which typically varies by a factor of 3-10 on the 
scale of Fermi energy. More recently, an anomalous be- 
havior has also been found in strong dependence of the 
low-temperature resistivity on the magnetic field parallel 
to the electron plane [[7[-]l4[. Since one expects such a 
field to couple only weakly to the orbital motion, this is 
usually interpreted as an indication that the electron spin 
plays an important role. Similarly as for the temperature 
dependence, the characteristic scale for the variations of 
resistivity with the parallel magnetic field seems again to 
be related to the Fermi energy ||ll| . 

Motivated by these intriguing experimental results, 
here I consider the 2D disordered Hubbard Hamiltonian 
as the simplest model of interacting disordered electrons, 
and propose a possible single explanation for both the 
temperature and the parallel magnetic field dependence 
of resistivity on the metallic side of the transition. Specif- 
ically, I assume that on the " metallic" side of the tran- 



sition the system is a disordered Fermi liquid, and study 
Boltzmann conductivity of the Hartree-Fock quasiparti- 
cles. First, I show that at zero magnetic field and at low 
temperatures screening of the random potential by the 
electron liquid enhances the Boltzmann conductivity by 
a factor of {1+g)'^, where g = UM/2 is the dimensionless 
interaction, and Af is the density of states at the Fermi 
level. Physically, this may be understood as that the ran- 
dom potential effectively became smoother, since at low- 
temperatures density of the electrons can adjust well to 
the given random configuration and screen it efficiently 
p^ . At high temperatures (and still at zero field), on 
the other hand, density is closer to uniform and more in- 
dependent of the randomness; consequently, the random 
potential is screened less, and the resistivity approaches 
it's (larger) value at g = 0. One point of this note is 
that the dependence on the parallel magnetic field can 
be similarly understood as a result of the interplay of the 
Hartree and the Fock terms in the screening process. It 
is primarily spin-down electrons that screen the random 
potential seen by the spin- up electrons (and vice versa), 
because the Pauli principle effectively reduces the inter- 
action between the particles with parallel spin through 
a competition between Hartree and Fock terms. This 
will be particularly true in the Hubbard model, where 
due to the assumed on-site repulsion the Hartree and the 
Fock terms for electrons with parallel spin actually cancel 
out. Crudely, at high fields (and low temperature) the 
electron system is strongly polarized and the number of 
down-spin electrons, and hence the screening, is reduced 
]l6[ . The second point is the expression for the Hall coef- 
ficient when the magnetic field is slightly out of the plane 
of electrons. In a recent experiment p^ the near absence 
of the field dependence of the Hall coefficient was used as 
an evidence against the present model: I show that this 
does not necessarily follow, and the two can be reconciled 
in a natural way assuming that the experimental sample 
is in the "clean" limit (see below). 

To make the above qualitative discussion more precise, 
consider the disordered Hubbard model for spin- 1/2 elec- 
trons on a quadratic 2D lattice: 

H = -t '^ c^^Cj^a +'^{vi + crH - /io)ni,o- + (1) 



(«J), 



U 






where c and c^ are the standard fermionic creation and 



annhilation operators. 



c' c is the number operator, 



and a = +,— labels the projection of spin along the 
direction of the magnetic field. I included the paral- 
lel magnetic field via Zeeman coupling to electron spin 
and set the value of the Bohr magneton and the Lande 
g-factor to unity, for simplicity. L'^ > represents the 
Hubbard on-site repulsion, and Vi is a Gaussian random 
potential with uiVj = Wdij, where the overbar denotes 
an average over randomness. To keep the algebra simple 
I will assmne low-filling, so that the electron dispersion 
is approximately quadratic, E{k) = k'^/2m, m = l/2i, 
and the Fermi surface is nearly spherical. The density 
of states per spin will accordingly be assumed to be con- 
stant M = m/2'K. Ho is the (bare) chemical potential. 

In the Hartree-Fock approach one searches for the opti- 
mal single-particle approximation to the above Hubbard 
Hamiltonian: 

Hhf = -t ^ c^^Cj^a +^ivi,cr + o-H - H)ni^^, (2) 

where the screened random potential is self-consistently 
determined as 



= Vi + 



U{r 



(3) 



and the thermal average in the last equation is taken 
over the states of Hhf- Chemical potential ^ is chosen 
so that Vi = 0. The screened random potential is local 
due to assumed on-site repulsion ||l8[. The last term in 
the Eq. 3 is the remaining difference between the Hartree 
and the Fock terms, which in case of the on-site repulsion 
and at -ff = differ simply by a factor of two. In gen- 
eral, for a given realization of the random potential the 
Hartree-Fock equations need to be solved numerically, 
which for a large system poses a somewhat non-trivial 
computational problem. Fortunately, for weak disorder 
the problem proves to be tractable analytically. To see 
this, expand the average on the right-hand side of the 
Eq. 3 to the first order in Vi. The Fourier components 
of the screened and the bare random potential are then 
readily found to be linearly related: 



Vaiq) 



v{q) 



i-c/n_,(g) 



0{v{q)' 



(4) 



where the (static) polarization function is given by the 
standard expression 



n,(g) = r^ 



dV 



1 



(27r)2 (iw„ - ^„{p)){iujn - S,a{p + 0)) ' 



(5) 



a.nd£,a{q) = E{q)-fi+aH. TheEqs. (4) and (5) describe 
screening of the weak random potential by the electron 
liquid. Neglecting the higher order terms in the Eq. (4), 
the screened random potential then satisfies 



Va{q}Va{p} = Waiq)S{q+p), 



(6) 




FIG. 1. Boltzmann resistivity (in units of Rg) as a function 
of parallel magnetic field (in units of Ef,o) at T/Ef,o ~ 0.1 
and for g — 2 for densities 0.7, 1, and 1.3 no. 
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Having obtained the two-point correlator Wa {q) for the 
screened random potential it becomes straightforward to 
calculate the Boltzmann dc conductivity of the HF quasi- 
particles [ |l9[ : 

where faiE) is the Fermi distribution function, and 

T-\E)^Af de{l - cos{e))W„{2V2mE sin -), (9) 

Ja 2 

is the scattering rate for the HF quasiparticles with en- 
ergy E. The chemical potential is fixed by the density of 
particles as n = J^J^a I fa-{E)dE. It is convenient to in- 
troduce a reference point: at a density no, the Fermi level 
is at Ef,q = fJ.{T = 0,H = 0) = no/2Af, and the Bohz- 
mann conductivity for g = is ctq — {e'^no)/(m2TrJ\fW). 
One can then write 



where, 



T 



2E 



F,0 V JO 



y — ^ — ia[y)dy, 
ay 



(1 - cos{9))de 



(10) 



" ii + Uo^f-^(y^Tsm\e/2W 

(11) 



Consider first few simple limits that may be treated an- 
alytically: 

1) H ^ Ep, T ^ Ep: one can neglect the ^-dependence 
under the integral in the denominator of Eq. (11), and 
the polarization function is approximately constant in the 
relevant region. It then follows that 
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FIG. 2. Fit to experimental data of ref. 11 for the den- 
sity 2.12 X 10^^ /cm^. Parallel magnetic field K is in units of 
the critical field (or Ef\ resistivity in 1000 Ohms, and the 
effective interaction is taken to be <; = 1.5. 



^ = (l+.9)2- + 0(cxp(-i?WT)). (12) 

(To n^ 

Conductivity is enhanced by the interactions since the 

random potential is effectively reduced by the electron 

liquid. 

2) H ^ Ep, T ^ Ep: the contribution of down-spin 

electrons to ub is exponentially small, as well as their 

screening of the random potential seen by the up-spin 

electrons. As a result 






no 



0(exp(-ff/T)). 



(13) 



3) if ^ Ep, T ^ Ep: the integral multiplying g in the 
Eq. 11 is of order Ep/T, so one can expand: 






no 



"(1 + 5 



n Ep^Q 
no T 



+ 0{{Ep.^o/Tf))). (14) 



The conductivity is a monotonic function of temperature, 
as a consequence of the assumed white-noise random po- 
tential (in contrast to ref. 15, where random Coulomb 
scatterers were assumed). 

At Fig 1. I display the low-temperature {T — Q.lEpfl) 
resistivity as a function of the parallel magnetic field. 
The saturation occurs around the corresponding Fermi 
energy, which is proportional to the density, and it ap- 
pears sharper than in the case of temperature depen- 
dence. Resistivity is also independent of the direction 
of the parallel field, and depends only on the degree of 
polarization of the electron liquid. The ratio between 
the saturation and the low-field value of resistivity is ap- 
proximately independent of density. All these features 
seem consistent with the experiments on Si MOSFETs 
[ pl| , iQ , iQ . The linear dependence of the resistivity 
on density in the plateau region, however, seems weaker 
than the observed fl^. This is also true at ii = 0, and 
presumably would be rectified by taking into account the 
long-range nature of Coulomb interaction |15| . 

On Fig. 2 a fit to data of Okamoto et. al. ||ll[ for 
the density 2.12 x 10^^/cto^ (critical density is at 0.98 x 



10^^/cm^, for comparison) and T = 0.21i^ is presented. 
The only fitting parameter is the effective interaction g 
which is determined by the ration of the plateau and 
the zero-field resistivities, and is taken to be g — 1.5 
(apart from the overall resistivity scale). To simplify the 
calculation I also set T = 0, which produced the kink 
at the critical field. The overall shape of the curve is, 
evidently, well reproduced. 

While the present calculation appears to describe the 
main qualitative feataures of the data, for a true quan- 
titative comparison with experiment it would need to 
be amended with the long-range Coulomb interaction, 
valley degeneracy, finite thickness of the electron layer, 
and weak-localization corrections. Since I assumed weak 
disorder, the present considerations are directly relevant 
only to the good metallic region [kpl ^ 1), where quan- 
tum corrections should indeed be negligible at not-too- 
low temperatures. Experimentally, however, even at the 
insulating side and close to the transition there is a simi- 
larly strong dependence of resistivity on the parallel mag- 
netic field [|2|. While this feature, strictly speaking, lies 
outside the reach of the present calculation, it would ap- 
pear less surprising if there was indeed no quantum phase 
transition at T = 0, since the system would then al- 
ways be in the same phase and would be expected to 
react to the external perturbation in a qualitatively sim- 
ilar way. 

After the completion of this work I learned of a pa- 
per by Dolgopolov and Gold where a similar mechanism 
for the magnetic field dependence of the resistivity was 
considered |jlq| . The principal difference with the present 
work (besides them assuming the Coulomb interaction) is 
in that Dolgopolov and Gold included only the Hartree 
term into screening, so that the spin-up particles pro- 
vided an additional constant contribution to the scatter- 
ing rates for both types of particles. Their conclusions 
are, nevertheless, in qualitative agreement with mine. 

While the obtained results seem broadly in agreement 
with the experiments on Si inversion layers, they describe 
the experiments on GaAs heterostructures less well [|3|. 
It is possible that the orbital effects of the type discussed 
in Ref. 20 are of greater importance there 122] . 

Finally, by assuming the magnetic field slightly off the 
electronic plane so that there is a weak perpendicular 
component i/j_ = Hsm9 <ti H, using the present two- 
species model one can straightforwardly calculate the 
Hall coefficient at low temperature: 



R 






A{H) 



mcB^{H) 
where the functions 



jIA^{H) 



A{H),B{H)=J2 



aBAH)ry(.H) 



(15) 



(16) 



and ll!j_ — eH±/nic is the corresponding cyclotron fre- 
quency, and (Tb^ct Boltzmann conductivities for particles 



with the spin projection a. For r_|_ = t_ the scatter- 
ing time cancels out, and one is left with the familiar 
Rxy/H±^ = 1/enc. In general case it does not, and 
the Hall coefficient in principle becomes field dependent. 
However, by expressing the Hall coefficient in terms of di- 
mensionless H/Hc, where He is the fully polarizing value 
of the field, one finds that the relative magnitude of the 
variation of the Hall coefficient with the field crucially de- 
pends on the numerical value of the parameter x = uJcTq, 
where cUc = eHc/mc and tq is the bare (unscreened) scat- 
tering time. For a; <C 1 one can neglect A{H) in the de- 
nominator of the Eq. (15) and the uj±Ta- factors in A{H) 
and B{H) to end up with pi) 



R 



xy 

Hi_ 



e a+{H)T+{H) + a^{H)T^{H) 
mc [a+{H) + a-{H)f 



(17) 



which is the expression used by Vitkalov et. al. ||I^ in 
the analysis of their measurements. In this case there 
should be a significant variation of the Hall coefficient 
with the field, unlike what has been found experimen- 
tally. This apparent contradiction between the theory 
and the experiment can be removed if one notices that 
in the opposite clean limit a; — > oo one can neglect B(H) 
in the Eq. 15, as well as unity compared to {uj^To-)'^ in 
A{H) to find that 



R 



1 



■xy 

H± nee ' 



(18) 



i. e. field independent once again. For an a; ^ 1 the 
Hall coefficient will thus be a weak function of the field 
(Fig. 3). A very crude estimate for the Vitkalov's exper- 
iments suggests a value of x of order ten for the density 
n — 1.44 X lO^^/crri^. For example, for a; = 30 and 
8 = 2.4° the variation of the T = Hall coefficient with 
the field is less that 5%. Finite temperature will addi- 
tionally decrease this number, so considering the overall 
level of approximation in this paper the theory and ex- 
periment appear to be in a reasonable agreement. 

In conclusion, I demonstrated that the Boltzmann con- 
ductivity of a disordered 2D Fermi liquid is in princi- 
ple dependent on the degree of it's polarization, and ar- 
gued that this may be a plausible explanation for the 
observed strong magnetoresistance of the electron liquid 
in Si-MOSFETs in a parallel magnetic field, in the low- 
temperature metallic phase. I calculated the Hall coeffi- 
cient based on this model and found it to be in qualitative 
agreement with a recent experiment. 
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by NSERC of Canada. 
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FIG. 3. The T = Hall coefficient normalized to its 7? = 
value as a function of the field H (in units of He) for = 2.4°. 
The upper curve is for the parameter x — and the lower for 
a; = 30 (see the text). The interaction is g = 1. 
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